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SUMMARY 


This  report  presents  a demonstration  of  the  usefulness  of 
the  ATLAS  system  in  performing  three-dimensional  elastic  stress 
analysis  of  a turbine  blade. 

Modeling  details  for  a shrouded  uncooled  turbine  blade  are 
outlined  and  program  execution  and  data  management  techniques 
are  discussed. 

It  was  concluded  that  three-dimensional  elastic  stress 
analysis  provides  an  accurate  means  of  predicting  stresses  in 
a complex  structure.  However,  high  computer  costs  require  that 
this  method  of  stress  analysis  be  used  vrith  discretion.  Areas 
for  further  study  are  suggested. 
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1.0  INTRDDUCTIOri 


Assessment  of  the  state  of  stress  in  solid  structures  has 
long  been  the  goal  of  stress  analysts.  Mov.'ever,  Lhr<^e-dimensionel 
stress  methods  have  been  limited/  for  the  most  pari , to 
photoelastic  methods  which  provide  good  results  for  a very  limited 
number  of  loading  types.  For  excimple,  cinalysis  of  turbine  blades 
has  been  limited  to  centrifugal  loading  because  thermal  and 
aerodynamic  loadings  have  been  difficult  or  impossil'le  to  sirulato 
oy  photoelastic  methods. 

In  recent  years  the  stiffness  method  of  finite  element 
stress  analysis  has  provided  a solution  for  three-dimensional 
stress  analysis  through  the  isoparametric  solid  element  (ref. 

1) , Since  1968  there  has  been  a proliferation  of  finite  element 
programs,  both  public  and  proprietary,  which  incorporate  some 
form  of  the  three-dimensional  isoparametric  solid  element.  The 
ATLAS  System  (ref.  2)  is  one  such  program,  available  to  government 
agencies  and  certain  of  their  contractors,  which  provides  the 
user  v/ith  a highly  versatile  isoparametric  brick  family. 

However,  users  soon  found  that  something  more  than  an 
accurate  finite  element  was  necessary  for  a successful  three- 
dimensional  stress  amalysis.  Even  relatively  simple  three- 
dimensional  models  can  produce  very  large  and  costly  comp\itational 
problems  which  may  exceed  the  capacity  of  the  largest  computers. 

It  became  evident  that  an  efficient  data  management  system  and 
substructuring  capability  were  as  inportant  as  the  finite  element 
itself  for  a successful  three-dimensional  stress  analysis. 

This  study  was  \indertaken  to  demonstrate  the  usefulness 
of  the  ATLAS  system  in  executing  an  elastic  stress  analysis  of 
a turbine  blade.  The  problem  selected  required  use  of  tlie 
system's  isoparametric  element  family,  various  loading  options, 
data  management  features,  and  automated  substructuring  capability, 
all  of  which  are  essential  for  successful  execution  of  i.ioderate 
to  large  three-dimensional  stress  analyses. 

The  authors  gratefully  acknowledge  the  support  of  the  program 
meinager,  M.  Aarnes  and  of  the  ATLAS  Staff  in  accomplishinn  the 
goals  of  tills  document. 


2.0  THE  ATLAS  SYSTEM 


ATLAS  is  an  integrated  structural  analysis  aiid  design  syster' 
operational  on  the  Control  Data  Corporation  (CDC)  6600/C'’IiER 
computers.  It  is  a modular  system  of  computer  codes  integrated 
within  a common  executive  cind  data  base  framework.  T!ie  systerr 
has  a broad  scope  in  that  its  analytica]  capabilities  support 
many  different  Put  related  aeroelastic  teclanological  disciplines . 
However,  in  this  report,  ATLAS  will  be  discussed  only  in  ..ts 
role  as  a three-dimensional  elastic  stress  analysis  tool. 

The  element  selected  for  this  demonstration  is  the 
isoparanveter  brick  element.  The  ATLAS  brick  element  famil”  is 
composed  of  four  major  elemeents  v/ith  orthotropic  material 
properties.  They  are  the  linear,  quadratic,  cubic  and  quartic 
bricks.  That  is  geometry,  displacement,  thermal  strains,  and 
pressures  are  expressed  as  linear  through  quartic  polynomials 
along  the  edges  of  the  brick.  Each  element  is  defined  by  3 
corner  nodes  cind  12  edges  each  of  vjhich  may  have  up  to  3 
intermediate  nodes.  Modes  are  allov’ed  3 translator^'  degrees 
of  freedom. 

Loadings  provided  by  ATLAS  include  any  or  al],  of  the 
foliating: 

a)  Inertia  loads 

b)  Point  loads 

c)  Pressure  loads 

d)  Thermal  loads 

e)  Specified  displacements 

Centrifugal  loading,  which  is  a special  case  of  the  inertia 
loading  option,  is  accomplished  by  defining  a rotation  vector 
which  provides  both  direction  and  angular  velocity  of  the 
rotation.  The  number  of  load  cases  in  a given  analysis  is  limited 
only  by  the  capacity  of  the  computer  system  used. 

Three  types  of  coordinate  systems  are  available  to  the  user. 
These  are  the  rectangular,  cylindrical  and  spherical  systems. 

Any  number  of  each  may  be  used  within  a problem  to  define  both 
input  and  analysis  reference  frames.  Thus,  complex  gecmotrics 
and  supported  and  specified  boundary  conditions  may  be 
conveniently  selected  by  the  user  in  order  to  simulate  nearly 
any  real  situation. 

Element  stresses  in  terms  of  3 normal  and  3 shear  stresses 
are  computed  for  the  global  coordinate  system  at  the  element 
centroid.  The  user  may  also  request  that  stress  at  tlx:  nodes 
be  printed  instead  of  or  in  addition  to  the  centroidal  stresses. 
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3.0  THE  BLADE  MODEL 


This  section  describes  the  geonetric  and  loading  models 
generated  for  a shrouded,  uncooled  pov/er  turbine  blade.  Blade 
geometry,  metal  temperature,  and  gas  pressures  were  provided 
by  Lycoming  Division  of  the  AVCO  Corporation.  The  subject  blade 
had  previously  undergone  a partial  three-dimensional  analys.is 
in  which  a portion  of  the  blade  V7as  modeled  three- dimensionally 
and  the  remainder  modeled  two- dine  ns  ion  a 1 ly . It  was  tlie  goal 
of  the  current  study  to  provide  a more  detailed,  fully  t)irae- 
dimensional  stress  analysis  of  the  blade  and  thereby  demonstrate 
the  usefulness  of  the  ATLAS  system  in  three-dimensional  stress 
analvsis . 


The  original  three-dimensional  model  of  the  i>lade  included 
that  portion  of  the  blade  from  .5  inch  above  the  platform  do^rn 
through  the  first  serration  of  the  fir  tree.  Torce  boundary 
conditions  from  a two-dimensional  model  >.;ere  applied  to  the  tO’-> 
of  the  irtodel  cind  support  boundary  conditions  v;ere  applied  to 
t!ie  fir  tree. 


'^he  current  model  was  built  from  the  original  by  adding 
the  remainder  of  the  fir  tree,  the  airfoil,  and  the  tip  shroud. 
Element  comer  node  numlDers  from  the  original  three-dimensional 
model  were  retained  in  the  updated  model,  but  the  element 
definition  of  the  original  section  \;as  upgraded  from  2D  or  24 
node  bricks  to  32  node  ’ -ricks.  The  purpose  of  the  refinement 
was  to  more  accurately  determine  the  stresses  in  the  area  of 
the  blade  root  fillet.  The  addition  to  the  fir  tree  was  modeled 
with  32  node  bricKs  and  the  airfoil  and  tip  shroud  were  nodolfjd 
with  8 node  bricks.  The  full  blade  model  v\?as  then  subs tructurcd 
as  shown  in  figure  3.1 


The  decision  to  use  substructures  \/as  made  at  tlie  cutset 
in  order  to  allov/  for  suljsequent  modifications  to  boundary 
conditions  and  tip  loads.  However,  it  was  found  during  -.-  recution 
of  tlie  problem  that  substructuring  would  have  Ijeen  necessary 
in  any  case  to  reduce  the  problem  size  to  fit  the  production 
configuration  of  the  Boeing  Computer  Services  (BCS) , CDC  66Q0 
computer.  That  problem  .is  discussed  in  more  detail  in  section 
4.0. 


Each  substructure  was  defined  by  a corresponding  stiffness 
data  set  containing  nodal  point  coordinates  and  element 
definitions.  The  substructure  numbers  and  data  set  numbers 
together  with  other  pertinent  data  ar^  given  in  table  3.1. 
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Figure  3.1  Turbine  Blade  Geometric  Model 


Table  3.1  Substructure  Oescripticn 


i 

o 


£ 'obstruct  ure  number 

11 

13 

14 

15 

16 

21 

Stiffness  data  set 

1 

3 

4 

c 

V 

6 

-- 

Number  of  nodes 

1076 

219 

57 

895 

953 

315 

.lumber  of  elements 

92 

81 

23 

66 

9 0 

-- 

’Nodes  per  element 

32 

8 

8 

32 

32 

-- 

Retained  nodes 

145 

46 

’6 

269 

154 

— 

Free  freedoms 

2616 

510 

120 

1869 

2 388 

9 4 5 

Retained  freedoms 

435 

138 

48 

807 

462 

-- 

Supported  freedoms 

132 

0 

0 

n 

-- 

Average  half  bandwidth 

848 

86 

76 

x281 

732 

531 

riaps  of  each  slabs tructure  together  with  element  corner  not]e 
numbers  are  given  in  figures  3.2  and  3. 4-3. 6.  Element  numbers 
in  data  set  1 were  generated  by  adding  the  tv/o  least  significant 
digits  from  one  of  tire  node  numbers  on  the  element  upper  face 
to  an  integer  equal  to  or  greater  than  4200.  This  nethor!  may 
be  demonstrated  by  observing  the  fir  tree  section  map  sh.own  in 
figure  3.3  in  conjunction  v/ith  figure  3.2.  One  may  obtain  any 
of  the  element  numbers  in  data  set  1 by  usi.ng  the  prefixes  given 
in  figure  3.2  and  adding  the  integers  found  encircled  in  figure 
3.3.  Element  numbers  for  data  sets  5 and  6 correspond  to  a node 
number  usually  found  in  the  second  quar.rant  of  the  element  upper 
face . 

The  stiffness  data  set  information  was  checked  for  accuracy 
via  the  plot  postprocessors  found  in  the  ATlu\S  system.  Additional 
checking  for  some  of  the  more  complex  geometries  was  done  on 
Vector  General  3-D  Vector  scope  driven  by  a PDF  11/45  ccn;iuter. 

The  global  coordinate  system  used  for  nodal  point  defirition 
was  a right-hand  rectangular  coordinate  system  witli  the  ;:-exis 
as  the  center  line  of  the  engine  (positive  aft) , the  y-axis  as 
the  tangential  direction,  and  the  z-axis  as  the  blade  stacking 
axis.  Two  additional  coordinate  systems  were  used  for  paiposcs 
of  boundary  condition  specifications.  The  x-axis  of  the  systems 
lay  parallel  to  the  longitudinal  axis  of  the  fir  tree  (positive 
aft).  The  pressure  side  system  was  rotated  45°  about  its  x-uxis, 
while  the  suction  side  system  was  rotated  -45°  about  ius  x-axis. 
These  systems  appear  in  the  ATLAS  coordinate  system  definitxons 
as  ROOTPS  and  ROOTSS  respectively,  and  are  shown  in  figure  J.8. 
Note  that  the  fir  tree  lands  were  free  to  slide  parallel  tc  the 
contact  surface.  IJo  shearing  forces  were  allowed. 

The  tip  shroud  boundary  condition  aJ lowed  translations  in 
all  directions  but  no  rotation  about  the  global  z-axis.  This 
was  accomplished  by  attaching  8 massless  beams  witii  stiffness 
properties  defined  in  the  x-y  plane  only  from  each  of  the  corners 
of  the  tip  shroud  to  2 nodal  points  on  the  upper  and  lower 
surfaces  of  the  tip  shroud  at  the  intersection  with  the  z-axis 
as  shown  in  figure  3.9.  The  2 nodal  points  were  constrained 
against  rotation  about  the  z-axis. 


REV  SYM 


NO.  D6-42735 

7 

Figure  3.2  Stiffness  Data  Set  1 


Figure  3.6  Stiffness  Data  Set  3 


Four  basic  loadcases  were  considered  in  this  analysis. 

They  were  centrifugal  loads,  thermal  loads,  aerodynamic  loads, 
and  a tip  rxob  load.  In  addition  to  the  basic  loadinas,  various 
combinations  of  the  basic  loads  were  considered. 

Centrifugal  loads  were  generated  by  an  external  data 
generator  since  the  ATLAS  inertia  loading  capability  v/as  not 
operational  at  the  time  of  this  analysis.  The  loads  generated 
were  due  to  rotation  at  16,000  rpm  around  the  global  x-axis  and 
were  input  to  the  ATLAS  system  as  external  nodal  forces. 

The  metal  temperatures  provided  by  Lycoming  are  shov/n  in 
table  3.2.  All  temperatures  input  to  the  program  were 
interpolated  from  the  Lycoming  data  by  a computer  program  in 
order  to  obtain  the  necessary  accuracy  for  the  fine  element  mesh. 


Table  3.2  Temperature  vs.  Z-Coordinate 

Temperature,  F Z-Coordinate , Indies 


1145 

4.00 

1145 

4.95 

1145 

5.00 

1150 

5.20 

1170 

5.70 

1192 

6.00 

1240 

6.50 

1285 

7.00 

1323 

7.50 

1354 

8.00 

1367 

8.50 

1362 

9.00 

1325 

9.50 

Aerodynamic  pressure  loads  were  also  provided  by  Lycoming 
and  are  given  in  teible  3.3.  Pressure  loads  were  computed  as 
a function  of  chord  length  and  Z-coordinate  and  input  to  the 
program  as  element  surface  pressure  loads. 

The  tip  rub  load  was  assumed  to  be  50  pounds  acting  in  the 
negative  tangential  direction  at  the  tip  shroud. 
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Table  3.3  Blade  Pressures 


Z=5.13  Inches 


Percent 

Pressure , 

psi 

Chord 

Pressure  Side 

Suction  S 

0 

17.87 

17.87 

10 

15.11 

12.54 

20 

15.20 

12.03 

30 

15.26 

11.95 

40 

15.26 

12.02 

50 

15.26 

12  .09 

60 

15.01 

11.99 

70 

14.90 

11.9  7 

80 

14.  79 

12.18 

90 

14.59 

12.78 

10  0 

17.87 

17.87 

Table  3.3  Blade  Pressures  (Continued) 


Z=7.30  Inches 


Percent 

Pressure , 

psi 

Chord 

Pressure  Side 

Suction  Side 

0 

19.90 

19.90 

10 

18.51 

15.37 

20 

18.44 

14.75 

30 

18.32 

13.95 

40 

17.95 

12.68 

50 

17.71 

11.83 

60 

17.61 

11.52 

70 

17.41 

11.75 

80 

17.01 

12.37 

90 

16. 33 

13.53 

100 

19.90 

19.90 

Table  3.3  Blade  Pressures 

(Concluded) 

Percent 

Chord 

Z=9.04  Inches 
Pressure,  psi 

Pressure  Side  Suction  Side 

0 

23.28 

23.28 

10 

19.57 

14.31 

20 

19.19 

13.66 

30 

19.10 

12 .65 

40 

19.15 

11.74 

50 

19.03 

11.  20 

60 

18.76 

10.99 

70 

18.18 

11.36 

80 

17.41 

11.70 

90 

16.82 

12.85 

100 

23.28 

23.28 
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4.0  PROGRAM  EXECUTION 


Execution  of  a stress  analysis  using  the  ATLAS  System  can 
be  as  simple  as  inputing  the  words  "PERFORM  STRESS."  ilouever, 
large  three-dimensional  analyses  will  require  more  user 
interaction.  For  those  analyses  requiring  such  interaction, 

ATLAS  provides  a concise  control  language  that  allows  step  by 
step  management  of  the  solution  process.  This  approach  u’as  used 
extensively  in  the  execution  of  the  demonstration  problem. 

The  solution  steps  for  the  demonstration  problem  v/ere  grouped 
into  five  logical  blocks  as  follows: 

1.  Input  data 

2.  Generate  stiffness  and  loads  matrices 

3.  Interact  substructures 

4.  Merge  and  reduce  substructures 

5.  Back  substitute  for  displacement  solutions,  calculate 
stresses  and  print  out  nodal  point  stresses. 

The  first  three  blocks  were  executed  for  all  substructures 
at  once  while  the  last  two  blocks  were  executed  once  for  each 
substructure  except  substructures  13  and  14  which  v;ere  executed 
together.  Each  of  the  above  blocks  utilized  one  or  more  ATLAS 
executive  statements. 

It  was  found  in  the  first  attempt  to  execute  the  problem 
that  the  initial  substructuring  arrangement  contained  a 
substructure  which  produced  more  data  than  could  be  stored  on 
a single  disk  storage  device  on  tlie  production  configuration 
of  the  DCS  CDC  6600.  The  job  would  have  aborted  due  to  a track 
limit  error.  The  offending  substructure  was  broken  into  two 
substructures,  numbered  15  eind  16,  and,  as  a result  of  this 
experience,  guidelines  were  set  up  which  reduced  tlie  probability 
of  further  track  limit  aborts.  These  guidelines  involved  careful 
data  management  through  the  use  of  substructuring,  ATLAS  executive 
statements,  cind  CDC  6600  job  control  cards,  together  v;itli  a good 
understanding  of  the  substructure's  gross  stiffness  matrix 
effective  half  bandwidth.  The  guidelines  are  outlined  lielow. 

First,  if  the  length  of  any  one  substructure  data  file  such 
as  the  stiffness  matrix  file  or  merge  data  file  exceeds  one-half 
the  capacity  of  the  device  on  which  it  is  stored,  the  problem 
size  should  be  reduced  by  replacing  that  substructure  with  two 
or  more  substructures  whose  data  files  fit  the  above  criterion. 


Second,  whenever  possible,  assign  the  largest  data  files 
to  different  disk  storage  devices  through  CDC  6600  job  control 
cards.  For  example,  restart  files  such  as  SAVESSF  should  be 
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assigned  to  a different  disk  than  the  files  being  loaded  to  or 
from  SAVESSF  such  as  STIFRNF,  etc. 

Third,  use  the  ATLAS  "SA’/E  flATRIX"  option  to  separate  data 
files  according  to  when  they  are  needed  in  the  execution  process. 
For  example,  store  the  element  stiffness  matrices  and  the  element 
stress  matrices  which  are  generated  at  the  same  time  on  separate 
save  files  so  that  they  may  be  loaded  as  needed  rather  t!ian  in 
one  large  file. 

Execution  of  tlie  demonstration  problem  without  substructuring 
would  have  required  up  to  40  million  v;ords  of  disk  storage  space. 
However,  with  the  use  of  the  guidelines  given  above,  the  required 
disk  space  was  reduced  to  about  12  million  words.  The  maximum 
lengtli  of  a single  file  was  about  3 million  words  which  v;as  well 
within  the  10  million  word  capacity  of  a single  storage  device. 
Thus,  the  problem  could  be  executed  on  the  production 
configuration  of  the  DCS  GDC  6600. 

The  interested  reader  is  referred  to  the  ATLAS  control 
progrcim  listings  in  appendix  A for  details  of  data  management. 
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5.0  STRESS  CONTOUR  PLOTS 


Contour  plots  of  the  pressure  and  suction  surface  radial 
stresses  were  made  for  stiffness  data  set  6.  Plots  for  each 
of  the  following  loadcases  are  shown  in  figures  5.1-5.12. 

a)  Centrifugal  force  (fig.  5.1,  5.2) 

b)  Aerodynamic  loads  (fig.  5.3,  5.4) 

c)  Thf  -mal  loads  (fig.  5.5,  5.6) 

d)  Tip  rub  loads  (fig.  5.7,  5.8) 

e)  Combination  of  cases  a-c  (fig.  5.9,  5.10) 

f)  Combination  of  cases  a-d  (fig.  5.11,  5.12) 

The  periodic  islands  that  appear  in  sonve  of  the  plots  are 
a result  of  the  contour  plotting  method  used.  The  magnitude 
of  the  stresses  at  the  center  of  those  islands  is  correct,  but 
the  adjacent  contours  are  distorted  as  far  as  location  is 
concerned  and  should  be  judged  accordingly. 
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Stress  contours,  ksi 


Stress  contours,  ksi 


Figure  5.12  Radial  Stress,  CF  + Aero  + Thermal  + Rub  Loads,  Suction  Side 
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6.0  TWO-  AND  THREE- DIMCNS ION AL  ANALYSIS  COMPARISONS 


The  subject  turbine  blade  had  undergone  several  earlier 
analyses  including  a two-dimensional  NASTRAN  plate  analysis, 
a three-dimensional  ATLAS  analysis  with  force  boundary  conditions 
provided  from  the  NASTRAN  analysis,  and  a full  scale  photoelastic 
analysis.  In  addition  to  the  analytical  data,  Lycoming  test 
experience  provided  information  as  to  where  the  blade  failed 
\jnder  overspeed  conditions.  These  data  are  compared  in  this 
section . 

Boundary  conditions  for  the  numerical  analyses  described 
above  varied  v/idely  according  to  the  complexity  of  the  analvsis. 
The  NASTRAN  plate  analysis  boundary^  conditions  assumed  that  the 
airfoil  was  rigidly  fixed  at  the  root  and  constrained  against 
rotation  at  the  tip  as  shov/n  in  figure  6.1.  The  combination 
analysis  using  NASTFAII  cind  ATLAS  assumed  that  the  fir  tree  was 
supported  against  radial  motion  at  the  bottom  of  the  first  land 
of  the  fir  tree  and  against  translations  and  rotations  along 
the  sides  of  the  fir  tree  as  shown  in  figure  6.2.  Note  that 
the  effects  of  the  upper  portion  of  the  airfoil  and  tip  constraint 
were  generated  by  NASTRAN  as  element  forces  and  subsequently 
input  to  ATLAS  as  nodal  forces. 

Boundary  conditions  for  the  photoelastic  study  were  provided 
by  a real  hardware  disk.  Since  the  disk  material  was  far  more 
rigid  than  the  plastic  blade  model,  the  boundary  conditions  for 
the  photoelastic  study  should  have  been  very  similar  to  those 
shown  in  figure  3.8. 

Field  experience  for  failures  due  to  overspeed  indicated 
that  the  most  frequent  point  of  failure  was  about  . 15  inch  above 
the  blade  platform.  This  location  corresponds  to  a Z-dimension 
of  5.13  inches  and  v;as  used  as  a basis  for  comparisons  of  the 
results  of  the  three  analyses  mentioned  above  and  the 
demonstration  problem  results.  Figures  6.3a  and  6.3b  show  the 
radial  stresses  due  to  centrifugal  forces  predicted  by  the  four 
methods  (at  a Z-dimension  of  5.13  inches)  for  the  pressure  and 
suction  surfaces  of  the  blade  respectively. 

The  NASTRAN  2-D  amalysis  predicts  average  stresses  well, 
but  it  does  not  follow  the  local  stress  shape  displayed  by  the 
photoelastic  model  well.  The  combination  3-D  and  2-D  analysis 
seems  to  have  similar  characteristics  to  the  2-D  run  only  greatly 
amplified  and  is  the  farthest  from  the  photoelastic  results. 

The  force  boundary  condition  is  thought  to  be  a major  contributor 
to  this  behavior.  The  fully  3-D  ATLAS  aunalysis  has  very  similar 
shape  characteristics  to  the  photoelastic  analysis,  although 
the  amplitudes  do  not  match  exactly. 

The  unknowns  in  the  photoelastic  study  probably  outweigh 
those  in  the  ATLAS  analysis.  For  excunple,  tJie  photoelastic 
material  has  a modulus  to  density  ratio  of  10^  while  the  71 3C 
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Tip  shroud  free  to  translate 
no  rotation  about  z 


Infinitely  riaid  platform  boundary 


Figure  6.1  Boundary  Conditions  for  NASTRAN  2-D  Analysis 


Boundary  forces  from  NASTRAN 


Figure  6.2  Boundary  Conditions  for  ATLAS  3-D  + NASTRAN  2-D  Analysis 
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blade  raaterial  has  a ratio  of 


Thus,  the  nonlinear  stifl 


effects  produced  by  the  centrifugal  forces  should  produce 
different  deflection  patterns  and  different  stresses  in  .ich 
model.  There  is  also  the  possibility  that  the  photoelastic  rode] 
u.iderwent  local  yielair.g  at  the  point  of  stress  concentration 
v/hich  would  change  the  stress  picture  considerably.  It  is 
interesting  to  note  that  if  yielding  had  taken  place  in  the 
region  of  60-80%  chord  on  the  pressure  side  of  the  pJiotoelastic 
model  the  stress  on  the  suction  side,  being  beiovJ  the  elastic- 
limit,  v.'ould  have  increased  due  to  the  transfer  of  load.  If 
one  applies  this  hypothesis  to  the  ATLAS  3-D  data  in  fig-ure  6.3, 
the  resulting  curve  shapes  would  be  nearly  identical  to  the 
photoelastic  curves. 

A comparison  of  the  ATLAS  3-D,  the  photoelastic,  and  test 
experience  data  is  shov/n  in  figure  6.4.  Assuming  that  the 
overspeed  failures  were  due  to  centrifugal  force  overload 
exclusively,  the  ATLAS  3-D  analysis  comes  much  closer  to 
predicting  the  failure  point  than  the  photoelastic  stutly. 

It  would  appear  from  the  results  of  t'.ese  comparisons  and 
other  related  experience  with  the  three-dimensional  ifsopa/anetric 
elements  that  they  provide  the  most  accurate  and  most  versatile 
three-dimensional  method  of  stress  analysis  available.  The 
analysis  of  any  geometry  should  be  feasible  through  the  use  of 
suLstructuring,  and  the  accuracy  of  the  results  shouJd  be  limited 
only  by  the  analyst's  ingenuity  and  experience  and  the  c;cmputcr 
resources  availcible. 
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7.0  CONCLUSIONS  AND  RECO'lMniJDATIONS 
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The  following  conclusions  are  a result  of  this  study: 

1.  ATLAS  provides  an  accurate  three-dimensional  elastic 
stress  analysis  capability  for  analyzing  complex 
structures  such  as  turbine  blades. 

2.  Efficient  management  of  very  large  volumes  of  data 
is  as  important  to  successful  3-D  analysis  as  the 
element  itself. 

3.  Substructuring  is  essential  to  three-dimensional 
analysis  of  complex  structures  such  as  turbine  blades 
and  disks. 

The  following  are  recommended  areas  for  further  study: 

1.  A nonlinear  analysis  using  an  iterative  displacement 
and/or  differential  stiffness  approach  should  be  made 
to  determine  the  effect  of  centrifugal  stiffening  on 
local  stresses. 

2.  Plastic  analyses  through  the  use  of  substructuring 
and  piecewise  linear  analysis  should  be  explored. 

3.  Known  stress  concentrators  such  as  fillets  should  bu 
modeled  with  different  order  bricks  and  results 
tabulated.  This  could  be  accomplished  by  modeling 
the  concentrator  configurations  documented  by  Roark 
(ref.  3)  or  Peterson  (ref.  4). 

4.  Boundary  condition  options  which  allow  tensile  or 
compressive  reactions  but  not  bot'n  should  be  tried 
in  modeling  fir  trees  in  order  to  study  the  effects 
of  loss  of  contact  along  the  lands. 
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APPENDIX  A — ATLAS  OUTPUT 


This  appendix  contains  a complete  listing  of  the  ATLAS 
execution  of  the  demonstration  problem.  It  is  intended  to  provide 
detailed  information  for  those  wishing  to  carry  on  further  studies 
of  tlie  demonstration  blade  and  to  provide  a general  guide  for 
anyone  intending  to  do  three-dimensional  stress  analysis  using 
the  ATLAS  system. 

The  microfiche  file  at  the  end  of  this  appendix  is  organized 
in  blocks  according  to  the  ATLAS  executive  modules  used,  and 
corresponds  to  the  following  outline: 

A.l  Read  input 

A. 2 Execution  of  STIFFNESS  and  LOADS  modules  and  printing 
of  stiffness  data 

A. 3 Execution  of  INTERACT  preprocessor  and  printing  of 
interact  data 

A. 4 Execution  of  SS-MERGE  and  SS- REDUCE  procedures  to 

perform  matrix  merge  and  reduce  operations,  and  SS- 
PARTITIOM  procedure  to  partition  the  highest  level 
substructure 

A. 5 Execution  of  SS-BACK  procedures  and  STRESS  and  OUTPUT 
modules  and  the  nodal  stress  output  routine 


Table  A.l 

Computation 

Time  and 

Resources  * 

Substructure  'lumber 

11 

13,14 

Procedure 

Time* , sec 

CRUs 

Time,  sec 

CRUs 

CF  Loads 

— 

54 

-- 

22 

Stiffness/Loads 

1696 

278 

89 

25 

tie  rge/ Reduce 

3673 

998 

35 

29 

Bac)csub/S  tress 

973 

254 

220 

6 3 

Totals 

6342 

15  84 

384 

139 

Table  A.l  Computation  Time 

and  Resources*  (Continued) 

Substructure  Mumher 

15 

16 

Procedure 

Time,  sec 

CRUs 

Time,  sec 

CRUs 

CF  Loads 

— 

46 

— 

46 

Stiffness/Loads 

1149 

19  2 

1551 

254 

Ilerge/'.leduce 

6667 

2026 

3034 

778 

BacksuL/S tress 

761 

214 

854 

199 

Totals 

3577 

2478 

5439 

1277 

Table  A.l  Computation  Time  and  Resources*  (Concluded) 

Substructure  ilumber 
21 

Procedure  Time,  sec  CRUs 

Merge/Reduce/Partition  802  156 

Problem  Totals  - 21554  seconds  and  5634  CRU's. 

*Dollar  costs  are  directly  proportional  to  computer  resource  units 

The  total  input  card  count  for  the  demonstration  proJ^l-'in 
was  approximately  10,900  cards. 
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A.l  Input 


A.l  Input  (Continued)  ^ 


A. 2 Stiffness  and  Loads  --  SS 
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A 2 Stiffness  and  Loads  --  SS  13,  14  Appendix  D6-4273j 
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A. 2 Stiffness  and  Loads 

A. 2 Stiffness  and  Loads  --  SS  16  Appendix  D6-4-i735 


A. 3 Interact 


A. 3 Interact' (Continued) 


A. 3 Interact  (Continued) 


A. 4 Merge  and  Reduce 


A. 4 Merge  and  Reduce 


A, 4 Merge  and  Reduce 


A. 4 Merge  and  Reduce 


A. 4 Merge  and  Reduce 
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APPENDIX  B — STRESS  CONTOUR  PLOTS 


Stress  contours,  ksi 


Figure  B.l  Equivalent  Stress,  CF  Load,  Pressure  Side,  Set  6 
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Figure  B.2  Equivalent  Stress,  CF  Load,  Suction  Side,  Set  F 
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Stress  contours,  ksi 


Figure  B. 
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Stress  contours,  ksi 


Figure  B.4  Equivalent  Stress,  Aero  Load,  Suction  Side,  Set  6 
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Figure  B.6  Equivalent  Stress,  Thermal  Load,  Suction  Side,  Set  6 
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Stress  contours,  ksi 


Squre  B.ll  Equivalent  Stress,  CF  + Aero  + Thermal  + Rub  Loads,  Pressute  Side,  Set  6 
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Figure  B.21  Equivalent  Stress,  CF  + Aero  + Thermal  Loads,  Pressure  Side,  Set  3 


D6-427aS 


KEV  SYM 


c/i 


DISTRIBUTION  LIST 


i 

i 


! 

i 


No.  of 

Copies  To 


Commander,  US  Army  Aviation  Systems  Command,  P.O.  Box  209, 

St.  Louis,  Missouri  63166 
10  ATTN:  DRSAV-EXT 

1 DRSAV-FE  (Cliff  Sims,  Maint  Engr) 

1 DRSAV-EQ  (C.  Crawford,  Sys  Dev  5 Qual) 

1 DRSAV-FES  (H.  Bull,  Corpus  Christ!) 

2 DRSAV-ZDR  (Ref  Library) 

Project  Manager,  Advanced  Attack  Helicopter,  P.O,  Box  209, 

St.  Louis,  Missouri  63166 
1 ATTN:  DRCPM-AAH,  TM 

Project  Manager,  Utility  Tactical  Transport  Aircraft  System, 

P.O.  Box  209,  St.  Louis,  Missouri  63166 
1 ATTN:  DRCPM-UA-T 

Project  Manager,  CH-47  Modernization,  P.O.  Box  209, 

St.  Louis,  Missouri  63166 
1 ATTN:  DRCPM-CH47M 

Project  Manager,  Advanced  Scout  Helicopter,  P.O.  Box  209, 

St.  Louis,  Missouri  63166 
1 ATTN:  DRSAV-SIA 

Project  Manager,  Aircraft  Survivability  Equipment,  P.O.  Box  209, 
St.  Louis,  Missouri  63166 
1 ATTN:  DRCPM-ASE-TM 

Project  Manager,  Cobra,  P.O.  Box  209,  St.  Louis,  Missouri  63166 
1 ATTN:  DRCPM-CO-T 

Project  Manager,  Iranian  Aircraft  Program,  P.O.  Box  209, 

St.  Louis,  Missouri  63166 
1 ATTN:  DRCPM-IAP-T 

Commander,  US  Army  Materiel  Development  and  Readiness  Command, 
5001  Eisenhower  Avenue,  Alexandria,  Virginia  22333 
4 ATTN:  DRCRD-TE 

1 DRCLDC,  Mr.  R.  Zentner 

Director,  Eustis  Directorate,  US  Army  Air  Mobility  R5D  Lab, 

Ft.  Eustis,  Virginia  23604 
1 ATTN:  SAVDL-EU-TAS 

1 SAVDL-EU-SS,  Mr.  J.  Robinson 

Director,  Ames  Directorate,  US  Army  Air  Mobility  R5D  Lab, 

Ames  Research  Center,  Moffett  Field,  California  94035 
1 ATTN : SAVDL-AM 


1 


No.  of 
Copies 


To 


f 


I 


I 


1 


1 

1 


1 

1 

1 


1 

1 


1 

1 

1 


1 

1 

2 

1 

2 

12 


2 


2 


Director,  Langley  Directorate,  US  Army  Air  Mobility  R§D  Lab, 
Mail  Stop  266,  Hampton,  Virginia  23365 
ATTN:  SAVDL-LA 


Director,  Lewis  Directorate,  US  Army  Air  Mobility  R^D  Lab, 
21000  Brook  Park  Rd,  Cleveland,  Ohio  44135 
ATTN:  SAVDL-LE 


U.  S.  Army  Industrial  Base  Engineering  Activity,  Rock  Island,  Illinois  61201 
ATTN:  DRXIB-MT,  Mr.  James  W.  Carstens,  Chief,  Manufacturing  Tech.  Division 


Air  Force  Materials  Laboratory,  Manufacturing  Technology  Division 
Wright-Patterson  Air  Force  Base,  Ohio  45433 
ATTN:  AFML/LTM 

AFML/LTN 
AFML/LTE 

Commander,  US  Army  Electronics  Command,  Ft.  Monmouth,  New  Jersey  07703 
ATTN:  DRSEL-RD-P 

DRSEL-GG-DD 

Commander,  US  Army  Missile  Command,  Redstone  Arsenal,  Alabama  35809 
ATTN:  DRSMI-IIE  " 

Technical  Library 
DRSMI-RSM,  Mr.  E.  J.  Wheel ahan 

Commander,  US  Army  Troop  Support  Command,  4300  Goodfellow  Blvd., 

St.  Louis,  Missour:  63120 

ATTN;  DRSTS-PLC 

Commander,  US  Army  Armament  Command,  Rock  Island,  Illinois  61201 
ATTN:  DRSAR-PPR-IW 

Technical  Library 

Commander,  US  Army  Tank -Automotive  Command,  Warren,  Michigan  48090 
ATTN:  DRDTA-RCM.l 

DRDTA,  Research  Library 

Commander,  Defense  Documentation  Center,  Cameron  Station,  Building  5, 
5010  Duke  Street,  Alexandria,  Virginia  22314 

Hughes  Helicopter,  Division  of  Summa  Corporation,  M/S  T-419,  Centinella 
Avenue  and  Teale  Street,  Culver  City,  California  90230 
ATTN:  Mr.  R.  E.  Moore,  Bldg.  314 

Sikorsky  Aircraft  Division,  United  Aircraft  Corporation, 

Stratford,  Connecticut  06497 

ATTN:  Mr.  Stan  Silverstein,  Section  Supv.,  Manufacturing  Technology 


i 

i 

! 





No.  of 
Copies 


To 


Bell  Helicopter  Company,  P.O.  Box  482,  Ft.  Worth,  Texas  76101 
2 ATTN;  Mr.  P.  Baumgartner,  Chief,  Manufacturing  Technology 

Kaman  Aerospace  Corporation,  Bloomfield,  Connecticut  06002 
2 ATTN:  Mr.  A.  S.  Falcone,  Chief  of  Materials  Engineering 

Boeing  Vertol  Company,  Box  16858,  Philadelphia,  Pennsylvania  19142 
2 ATTN;  R.  Pinckney,  Manufacturing  Technology 

Detroit  Diesel  Allison  Division,  General  Motors  Corporation, 

P.O.  Box  894,  Indianapolis,  Indiana  46206 
2 ATTN:  James  E.  Knott,  General  Manager 

General  Electric  Company,  10449  St.  Charles  Rock  Road, 

St.  Ann,  Missouri  63074 
2 ATTN;  Mr.  H,  Franzen 

AVCO-Lycoming  Corporation,  550  South  Main  Street,  Stratford, 
Connecticut  06497 

2 ATTN;  Mr.  V.  Strautman , Manager,  Process  Technology  Laboratory 

United  Technologies  Corporation,  Pratt  6 Whitney  Aircraft  Division, 
Manufacturing  Research  and  Development, 

East  Hartford,  Connecticut  06108 
2 ATTN:  Mr.  Ray  Traynor 


2 

1 

1 

1 

1 

1 

1 

I 


Director,  Army  Materials  and  Mechanics  Research  Center, 
Watertown,  Massachusetts  02172 
ATTN:  DRXMR-PL 

DRXMR-AG 
DRXMR-PR 
DRXMR-CT 
DRXMR-X 
DRXMR-AP 
DCXMR-M 
Author 


M </l  W i/t  C -A 

ii!S  S 

< <Q  5 c 

_i  Q c t: 
o uj  o»  « « ^ 
zP  ^ ^ 

3 _ w 3 


U CM  Z O- 

u r»»  •—  'O  6 

^ £D  C o 
41  CM  0£  « O 

2 ® ? • 01  ^ ‘ 
C£  >—  C CM 

«/t  < 41  « ^ I 
«/l  4^  3 »—  CO 

O 4-> 

4»0  « ^ 
c <y«  on  •*- 
<a  a «/n  < c 

^ ^ ‘ o 

u o on  «t »—  ■♦-* 

41  «q  ^ lO  O) 

£ «/t  -I  • C 

«A<  tf  u*»- 
^ « 2 k z 

e £ < *4101 


r-  S UJ  41  O • 

m OOC  ^ 4) 

*J  H-  O'*— 

i.  t.  on  . e 4-> 

« 41  * 

4-*  4->  o 4>  « 


§ -c 

c-S:s 

•*-0  3 «0 

£ o • c 2 ^ 

£41  41  z >> 

X 4i  4-*  4->  -O 
«J  > 41  W U O 

(A  4>  £ **-  41  4^ 

>*TJ  e 4J  T3  £-  *A 

vt  « <0  41  . 

f-  w 44  w a i- 

i/>  ^ O'  <4  O.  O 41 

< o -C  41  £ 
_j  41  c *; 

c CL  o a 

«i  •*-  "O  o>  *•- 

^ T9  41  Ut  «•* 

41  u c *o  c Ot  u 

£ a <C  3 <4  o o 

^ 41  U 

u E 

«*-  <4  41  C i.  «/» 

o C O 41  41  <4 

U 40  41 

o 5 2 i: 

4/1  ^ 4lt  ^ Ok<4 

41  «n  a n>  a E 

c *r-  o 2 u o 

.—  i/t  u u * 


41  4Q  10  c c 

4rt  C <4 

a <4  41  • -CO) 

-D  T3  </»  t 

41  M <4  41  41  * LJ 

Z ot  ^ «n  O L-  trt 

40  4i  ^ 4/t  41  **- 

V.  a > > “O 

4^1  41  U O 41 
o »n  c o»  u S z 

Q.  O ■*0 

0 0.0*0  X 

O •*•  I-  ov  2 

•^  4J  3 4l  _ 

40  l/t  40  k.  0»  • *0 

(Q  «0  «Q  >l  «A  41 
t_  »—  O *—  41  o> 

40  4i  41  (D  14  k 3 

«/l  ^ 41  C 3 

c ^ o a 14  «o  o 

0 <4  o O U -A 

E c o -•-  cA  a 

S o c C «A  i-  «A 

•o  •*'“  ax:  41  *o  ••- 

(/)  u (-  «A  «A 

«4  C *0  4l  V>  >\ 

4l  4l  4^  «A  X /— 

01  E T)  4)  «4 

40  .X  a o ^ c 

C "O  O C Q.  <4 

41  • ^ ^ e 

ot  41  z E «/)  o «/t 

41  4l  4/1  41  ns  O «A 

. C 1.  O'*-  41 

O.Z  14  14  41  <4  i- 
4J  C 40 

40  1.  44  *-  c 01 

' t.  O O E « 

I o c V c «f-  *o 

; a.**-  <4  O lA  O 41 

41  E «A  4^  41  «•* 

fc.C^«iAOVTSt/> 

o T3  C (A  O 41 


32  i-S 

M OC  «A  O*  *4 

u. *0  *i  r: 

*-•(/)  k.  o>  c .0 

ss  $ 

Sc  <4  41  C 

_J  Q X C r—  ■** 

0 U1  4l  10  4l  ^ 
Z ►“  V; 

3 m id  41  a 


lA  4.  14  ft!  ^ S 
.^040*eE^'*^<^ 
X4i4iC'^«^ft!3 
h-  Q.TS  «4  T>  CA  6 0* 


*-  O x: 

<u  •*-  **  V 

C "O  40  t 

o a 0 *0  « 

£ o • c 14 

C 41  4i  •*-  £ X 

S X 01  40  40  t? 

40  • 41  ^ o a 

(A  41  x:  ••-  41  40 

XT5  6 *'*  *0  t iA 

1/1  14  <4  41  — 

*-J.40|.3i- 
I/S  .0  O <4  O-  O 41 

< ox:  u X 

_j  4i  k-  ♦0  C 40 

c a o a 

< .»-  -O  ‘AHi* 

.O  ^ Ql  lA  40 

4lUC^C<Ak- 
£ a <4  3 <4  o o 

4^  40  *—  4|  U 

■0  u E 

14  01  C W i/l 

0 c o 01  oi  « 

1*.  •*-  o 4-1  40  41 

i/>  O ■—  IQ  3 U 

iO  40  lA  k.  0.« 

01  vs  a <4  a £ 

c — o 2 u o 

M (A  ou* 

a X 41  «o  10  c 

It-  M k.  >—  x:  o 
^ 10  «e  c c ■*' 

lA  C 10  — 40 

a <4  41  • .C  41 

^ Tj  irt  1. 

41  tA  14  4l  Ol  ' 0 

f vs  ^ ^ k.  w> 

40  41  £ lA  4l  — 

u a > > -o 

<4.  40  41  O O 4l 

O vs  c vs  k-  2 Z 

o.  O 40 

c u ^ ^ * 'C 

O — fc-  vs  2 

M 40  a « ^ 

le  <4  <4  X vs  41 

k-  T3  »“  ft*  ^ 

40  O;  4i  lA  10  k.  a 

VS  ^ 4)  c a 

c — o a « «o  41 
o 10  o O’  o .0 

• E ff  u — vs  a 

S O C C lA  i-  vs 

I ^ -t-  a £ 41  *o 

vs  u k.  vs  lA 

! 14  C *D  01  40  X 

^ O 40  vs  X ^ 

» vs  E V Ol  <4 

» 40  t-  a 40  kl »—  c 

1 c "O  O C -*-  Q.  » 

.40  «»  • *“  i 40  g 

tk.  vsoiocEvsOvs 

>0  oiOivsSiokivs 

> Ck  u u o—  o> 

Ol  0.£i4i44in>k. 

<X  40  C 40 

“ 40  k.  14  *-  C VS 

Sf~  ue7-oE«  — . _ 

s ns  o c c ^ 

• c CL-r-  lOOlAOJ 

• •*.  4>EVS40.,-4*  4 

-It  k-Cr-i4VSVS^ 

o — t3  c VS  O 

vs  «4-  «0  __  d[  4>  £ 

•F*  w *0  t»  B C 40  I 
£ 0>  4>  C — 40  , 

W Q.T5  *4  *0  vs  6 ’ 


-h 


M QC  «A  V*  I4 

u.  ►-  X)  •*-  *0  •- 

M </s  k.  VS  c £ 

^;s  s « 

< 14  01  c 

«l  Q X C r-  — 

kJ  Ul  01  <4  Ol  £ 

XV-  X W 

3 M lA  41  a 


O CM  Z & 

*“  ►**  K 

14  00  C 0 

41  CM  oe  <4  kJ 

vs  o ? 

41  P • 41  4» 
OC  *—  C <M 

vs  «t  41  14  — 
vs  *0  a r“  9? 
u 40  it  E <xos 
1-  41 0 lO  *• 

C vs  I/S  — 

14  a vn  «t  c 
£ £ — • O 

k>  u vn  < *“  40 

is:; 

T,SS“ 

C £ «C  - 41 

” .^4:|a 

V*  c i/n  2 O 


40  ►-  o»*- 

U k.  S/s  • C 40 

41  41  _ £ — 40 

£ 40  O 41  10 


^ 40  41 

C *«  40  k. 

-r*  O a 0-40  .0 

£ k»  • c 14 

g 41  41  — £ X 

41  K Ol  40  40  -o 

40  • 41  k.  ki  a 

vs  41  £ •*-  01  40 

X 'O  £ *0  T3  k V> 

VS  14  <4  41-*- 

L.  40  i.  3 k. 
</S  £ O'  14  O.  O 01 

^ O £ « x: 

_J  01  W 40  I*-  k 40 

c a.  O a 

.f-  T3  SA14- 

£ T5  01  vs  40 

41  k.  C ^ C vs  k. 

£ a 14  a 14  o o 

40  40  ^ 0!  kl  s*- 

-o  u 6 

S«-  10  01  C k lA 

0 C o 41  Ol  14 

S^  •*>  U 40  40  41 

vs  o *-  <4  a k 

lA  40  vs  k 

01  vs  a s a E 

c — o 2 kl  o 
F-  sA  ki  <0  • 

a X 01  40  <0  c 

«*-  f—  k £ O 

41  14  10  C C — 

vs  C to  ••-  40 

3 «0  41  • £01 

•S3  "O  V*  k 

41  vs  <0  41  41  * kl 

£ VS  *—  vs  ^ k sA 

*0  41  £ VS  .*-  41  — 

k a > > V 

s»-  40  01  u O ft! 

O sA  r «A  k 2 £ 

a o 40 

c 0 £ XI  X 

O — k sA  2 

-*-  40  a 41  — __ 

40  vs  40  k vs  • 'O 

10  <4  0 X vs  41 

k.  — ■ XI  •—  41  vs 

40  4/  41  V*  <4  k a 

vs  r-  ft!  c a 

c *-  o a <4  40  ft! 

J<4  O O U £ 

c u — VS  a 

O C c lA  k sA 

ax  OP  40  -F- 

vs  u k vs  V* 

14  C *0  Ol  40  X 

gOl  *1  lA  X ^ 

^ Ol  « 

40  •*-  a 40  k»  — c 

c *0  O C •*-  Q.  « 
40  41  • k ft!  40  E 

k V>Ol£0VsO«A 

O 4iOiVSW0klvs 


40  I.  «o  .—  C «A 

kOOEio—  • 
jiOOCk  C k-X? 
-C  O'-  iwovsooi 


■F-  w 40  47  g k 40 

£ Ol  O'  C - 40  4| 
k-  Q TJ  ns  O sA  E 


vn  •—  O >>  ^ 


£ Ui  u.  14 

U CM  Z O- 

k-  ■-«  X>  E 

14  *-  A c o 

Ol  CM  Z « kO 

vs  O p _ 

41  P • 01  4» 

CX  C CM 

VS  < Ol  14  — 

lA  40  S'-® 

i>  40  u.  E o ^ 

- 01  o « k 
c vs  vn  •- 
14  a c/>  4 c 

£ £ •—  • o 

sj  kl  «A  ^ — 40 

Ol  <4  X « OS 
£ VS  _l  • — C 
VS  « Z O — 

? £ «*:  “ t ^ 

vs  c vn  X o 

M 2 UJ  01  k.1  * 

14  O Z -1  4l 
.F.  4OK  CT  — 
V.  w vn  • c 40 

01  01  X 40 
40  40  o fti  n 

2 5 A o m s^ 


c P 40 

— o a p 

£ ki  • I 

E Q>  01  1 

S K 41  4 

«0  • 41  k I 


k-  C O.  I 
«I  T-  X) 
£ *0  01  P 
01  k C XI  ! 

£ a <4  a ' 

40  40  Ml 
t3  u I 
«4-  14  0;  c 
O CO 


01  *4  >0  I 

VS  C « 

a 14  OP  ■ 
*0  *o  I 

01  vs  m 01  < 

£ vs  M vs  n 

40  01  £ vs  1 


•*-  40  3 O'  • 
40  vs  40  k. 
m 14  10 

k M XJ  * 
40  01  01  VS 
VS  *—  ft! 


« c X)  01  4 

ft>  01  4V  ( 

vs  E -O 

40  a 40  I 

c ^ o c 

41  1 k « 4 

vs  41  £ E • 

E4>  «A  OP  ■ 

k O'* 

CX  £ 10  14  ' 

40  c 

♦0  k 14  * 

k kJ  E ' 

o c •* 


V*  k 0 
- k 4J  T3 
£ O'  O'  C 
k CX  TS  14  ' 


•-«(/>  k.  «/t  c ^ 

(/>  M O >,  4J 

O » •-  B 41 

^ <0  4>  C 


.•Mkaw 
r « o.  O’  V 
) £ 4J  £ 
. 4-»  k-  C W 


3 0>  «/t  k> 

: o c «/i  k 

B 3 « o o 

r-  4»  U V 


<0  3 W 
■>  4rt  W 0.4 
a « 3 e 
3 2 O O 

(J  <J  > 

“4=  6 

“ S5i; 

j 4 f ^ 

3 0 «/i  k 

0 4)  0)  ‘ O 

- 4/)  O k t/i 
3 *rt  4>  •♦“ 

3 > > *0 
U U O ft 
: k I £ 

- 0 o -*-* 


^ O kO  « 

U.  k-  T> 

*-■  CO  k lA  C £ 

irt  — O K ft 

S°  * sf  g 

_l  Q »N  C r- 

O UJ  ft  « ft  £ 

*c  k 


< k 3 

£ •»  ft 

: ~0  k *A 

g 3 ^ 

r ft  o o ft 

> £ ft  £ 

. «-*  V-  k ^ 


3 ft  «/*  «>• 

: T3  C >A  k 

B 3 ft  O O 
^ ft  kl  4*- 
3 U E 
DC  k <A 
; O ft  ft  ft 
• u £ ft 
> ft  3 k 
i lA  k O 4 

55  36 

o u • 


j tA  k 
> ft  ft  * u 
- lA  -9  k lA 
3 4A  v>  ft 
3 > > O 


J k 4A  . -O 
ft  >,  «A  ft 
3 ^ ft  «/• 

B 4/1  ft  k 3 

- ft  C 3 

3 3 A k»  ft 
3 0 U £ 
j •«-  lA  3 
::  c <A  k 4/1 
a £ ft  kl  -r- 

U k 4A  (A 

3 ft  4->  >t 
B «A  X ^ 

3 ft  ft 

3 4^  u ^ e 

-I"!'’ 

: B «/>  o cA 

A ft  ft  U kA 
ft 

B ft  ft  ft  k 
C 4> 

k ft  r-  C «A 

D e ft*^ 

- C V-  tJ 
ft  O V)  O ft 

A 44  •»-  ft  44 

- ft  <A  V>  "O  4A 
> *0  C «A  O ft 
Q ft  ft  £ O 

4 “t3  e k 4.»  O' 


A >S  lA  ft 
! 4A  ^ 3 

■ ft  C 3 

> 3 ft  44  ft 

> O’  U £ 


' ft  lA  lA  *0  4A 
' T3  C kA  O ft 
} ft  ft  £ O 

* T3  6 k *4  o 

> C •*-  44  ft  3 

> ft  T3  kA  6 kA 


M a;  4A  4A  ft 

U.  ■53  •«-■  *4  •— 

in  k 4A  c £ 

k/3  — £ >k  ft 

k/i  O 3 ^ £ ft 

< ft  ft  c 

_J  O >,  C r- 

k)  UJ  ft  ft  ft  £ 

X *-  X k 

3 4A  ft  3 

X <A  44  44 

^ ft  ft 

^ k c <A 


» ft  ••-  £ >v 
; ft  4^  44  ^ 

I k O 3 


> 44  V.  3 k 
7 ft  O O ft 
3 £ ft  £ 

to  44  k»-  k 44 

3 0 3 

■O  kA  »•- 

3 ft  4A  44 
: T3  C »A  k 

B 3 ft  O O 
ft  kj  «4- 

3 kj  £ 

U C k «A 

o ft  ft  ft 

• ft  44  44  ft 
- ft  3 k 

4 «A  k Q.C 

3 A 3 E 
3 S O O 

ft  ft  . 


O 0>  ft  ft 

k4j  •—  E O 

•—  OC  kA  kA  A 

u.  ►-  O •—  *4  r- 

in  k kA  c £ 

— o >,  ft 

•no  X r-  E ft 

< ft  ft  c 

_j  o >.  c 

k4  U4  ft  ft  ft  £ 

ac  V-  X k 

:3  — lA  ft  3 


ft  ••>  44 

D • £ ft 

3 TJ  kA  k 

B ft  ft  • O 

• «A  O k »A 


> k kA  . ■£ 

ft  >.  kA  ft 

} P—  ft  kA 

> 4A  A k 3 

• ft  k 3 

> 3 A «4  ft 

I CT  £ 


B «4  lA  X ^ 
3 ft  A 

3 «j  k.'  ^ C 


^ ft  A U «A 

O’*—  ft 

B A ft  A k 

c 44 

- A A-  C «A 

3 E « -•' 

C V *0 

ft  O «A  O ft 

■»  44  ft  44 

“ ft  *A  O *A 

> O C 4A  O ft 

a ft  ft  £ cr 

^ E k 44  cr 

4 C - A 3 

) ft  o »A  e *A 


